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ABSTRACT 

In this paper, we present the first results from the Renaissance Simulations , a suite of extremely high- 
resolution and physics-rich AMR calculations of high redshift galaxy formation performed on the Blue Waters 
supercomputer. These simulations contain hundreds of well-resolved galaxies at z ~ 25 - 8, and make several 
novel, testable predictions. Most critically, we show that the ultraviolet luminosity function of our simulated 
galaxies is consistent with observations of high-z galaxy populations at the bright end of the luminosity func¬ 
tion (M 1600 < “17), but at lower luminosities is essentially flat rather than rising steeply, as has been inferred 
by Schechter function fits to high-z observations, and has a clearly-defined lower limit in UV luminosity. This 
behavior of the luminosity function is due to two factors: (i) the strong dependence of the star formation rate on 
halo virial mass in our simulated galaxy population, with lower-mass halos having systematically lower star for¬ 
mation rates and thus lower UV luminosities; and (ii) the fact that halos with virial masses below ~ 2 x 10 8 M 0 
do not universally contain stars, with the fraction of halos containing stars dropping to zero at ~ 7 x 10 6 M 0 . 

Finally, we show that the brightest of our simulated galaxies may be visible to current and future ultra-deep 
space-based surveys, particularly if lensed regions are chosen for observation. 

Subject headings: galaxies: formation - galaxies: high-redshift - galaxies: luminosity function 


1. INTRODUCTION 

The most fundamental characteristics of the earliest galax¬ 
ies are challenging to determine directly. These galaxies lie 
at the edge of observability, or beyond, for even the largest 
current ground- and space-based telescopes, and inferences 
based on local stellar populations (from, e.g, the Milky Way’s 
stellar halo) are uncertain at best. In particular, it is difficult 
to determine the characteristics of the luminosity function of 
high redshift galaxies, which directly influences how reion¬ 
ization proceeds in the early universe. A recent example of 
this are efforts to use the Hubble Ultra Deep Field and its ex¬ 
tensions (Beckwith et al. 2006; Finkelstein et al. 2014; Koeke- 
moer et al. 2013; Bouwens et al. 2015), and, separately, lensed 
observations of high redshift galaxies from the Hubble Fron¬ 
tier Fields (Atek et al. 2015; Ishigaki et al. 2015), to estimate 
the luminosity function for galaxies at z — 6-10. These obser¬ 
vations succeeded in directly measuring the stellar luminosity 
function for galaxies with luminosities within a few magni¬ 
tudes of L*, but no deeper. 

The work of Bouwens et al. (2015) and Atek et al. (2015) 
provides only weak constraints on the faint-end slope of the 
luminosity function, which has profound theoretical implica¬ 
tions (e.g., Robertson et al. 2010; Fontanot et al. 2014). If the 
slope of the luminosity function is very steep, this implies a 
vast number of dim galaxies, unobservable by current instru¬ 
ments, that can produce more than enough ionizing photons to 
complete and sustain reionization with no additional sources 
(such as high-redshift black hole populations) needed. If, on 
the other hand, the slope of the high-redshift luminosity func¬ 
tion is shallow, there will be fewer galaxies than one might 
expect, with too few ionizing photons and thus the need for 
additional sources of ionizing radiation. 


A second question related to the shape of the high-redshift 
luminosity function is that of its ending point: how bright are 
the smallest high-redshift galaxies? Simulations of Popula¬ 
tion III stars and the transition to metal-enriched star forma¬ 
tion suggest that the smallest halos to form stars have masses 
of around 10 6 - 10 7 M 0 (O’Shea & Norman 2007; Johnson 
et al. 2008; Smith et al. 2009; Wise et al. 2012b; Crosby 
et al. 2013), with correspondingly few stars and low lumi¬ 
nosities. These simulations’ predictive capabilities, however, 
suffer from challenges relating to small-number and small- 
volume statistics. 

An additional issue relating to the study of high redshift 
galaxies relates to cosmic variance. The James Webb Space 
Telescope will certainly be able to see substantial numbers of 
galaxies at z > 10, but has a small field of view (Gardner et al. 
2006). As a result, interpretation of any JWST survey must 
by necessity take into account cosmic variance (e.g., Trenti & 
Stiavelli 2008; Moster et al. 201 1). While this can be done by 
using many small, widely spaced fields instead of one larger 
field (Sibthorpe et al. 2013), this may be undesirable, and thus 
a deeper understanding based on theory may be necessary to 
correctly interpret high-z survey results. 

In this paper, we address these pressing issues regarding 
the high redshift galaxy luminosity function using the Renais¬ 
sance Simulations - a trio of physics-rich simulations of high 
redshift galaxy formation that resolve several hundred galax¬ 
ies apiece, and which explore a range of cosmic environments 
at z > 8. Using these simulations (described in Section 2), 
we show results (in Section 3) that predict the high redshift 
galaxy luminosity function will flatten at magnitudes that can 
be probed by the next generation of telescopes. We discuss 
the implications of this work in Section 4. 
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2. THE RENAISSANCE SIMULATIONS 

The simulations were carried out using Enzo 1 (Bryan et al. 
2014), an open-source adaptive mesh refinement code that has 
been extensively used for simulating cosmological structures, 
and in particular high-redshift structure formation (e.g., Abel 
et al. 2002; O’Shea & Norman 2007; Turk et al. 2009; Wise 
et al. 2012a,b; Xu et al. 2013; Wise et al. 2014). Notably, H 2 - 
photodissociating and ionizing radiation from stellar popula¬ 
tions is followed using the MORAY radiation transport solver 
(Wise & Abel 2011). The properties of hydrogen and helium 
are calculated using a 9-species primordial non-equilibrium 
chemistry and cooling network (Abel et al. 1997), supple¬ 
mented by metal-dependent cooling tables (Smith et al. 2009). 
Prescriptions for Population III and metal-enriched star for¬ 
mation and feedback are employed, using the same density 
and metallicity criteria as Wise et al. (2014) but with a Popu¬ 
lation III characteristic mass of 40 M 0 . 

We simulated a region of the universe 28.4 Mpc/h on a side 
using the WMAP7 best-fit cosmology. Initial conditions were 
generated at z = 99 using MUSIC (Hahn & Abel 2011), and 
a low-resolution simulation was run to z = 6 to find regions 
suitable for re-simulation. The volume was then smoothed 
on a physical scale of 5 comoving Mpc, and regions of high 
((S) = (p)/(Q M pc)~l — 0.68), average ((S) cx 0.09), and low 
((S) ~ -0.26) mean density were chosen for re-simulation. 
These subvolumes, designated the “Rare peak“Normal” 
and “Void” regions, with comoving volumes of 133.6, 220.5, 
and 220.5 Mpc 3 , were resimulated with an effective initial res¬ 
olution of 4096 3 grid cells and particles in the region of inter¬ 
est, giving a dark matter mass resolution of 2.9 x 10 4 Mq. 
We allowed further refinement based on baryon or dark mat¬ 
ter overdensity for up to 12 total levels of refinement (i.e., a 
maximum comoving resolution of 19 pc). For more simula¬ 
tion details, see Xu et al. (2013, 2014); Chen et al. (2014). 
These simulations were evolved to z = (15,12.5,and 8) for 
the Rare peak , Normal , and Void simulations. In these sim¬ 
ulations, the halo mass function is well-resolved down to 
~2x 10 6 M 0 (~ 70 particles/halo), and the simulations con¬ 
tained (822,758,458) galaxies having at least 1,000 particles 
(M vir ^ 2.9 x 10 7 M 0 ). The simulations were analyzed using 
the YT analysis tool (Turk et al. 201 1). 

3. RESULTS 

Figure 1 shows the ultraviolet (UV) luminosity function for 
the galaxies within the refined region of our simulations. This 
luminosity function is calculated by extracting the mass, for¬ 
mation time, and metallicity of all star particles in each simu¬ 
lated galaxy at the final timestep of each simulation, and using 
the stellar population synthesis models of Bruzual & Chariot 
(2003) to create a spectrum for that galaxy (with each star 
particle representing a star cluster with a Chabrier mass func¬ 
tion). The UV luminosity in a window of A A = 100 A cen¬ 
tered on A = 1600 A is extracted from the resulting spectra. In 
each panel, the redshift-dependent fit to the Schechter func¬ 
tion from Bouwens et al. (2015) is shown at the same redshift 
as the simulation data, with the shaded region representing the 
1 -<j uncertainty in their fit parameters. The figure shows the 
luminosity function derived for all three simulations at z = 15, 
the Normal and Void simulations at z = 12, and the redshift 
evolution of the Normal simulation from z = 18 to z = 12 (212 
to 375 Myr). In each panel, we have calculated the magnitude 

http://enzo-project.org/ 





Figure 1 . UV luminosity function for galaxies in the Renaissance Simula¬ 
tions for the Rare peak , Normal, and Void simulations. The x-axis shows the 
absolute UV magnitude (Migoo) and the y-axis shows galaxies per magni¬ 
tude per comoving Mpc 3 . Top row: all simulations at z = 15. Middle row: 
Normal and Void simulations at z = 12. Bottom row: Normal simulation at 
z = 18,16,14, and 12. In the top two panels, the redshift-dependent fit to the 
Schechter function from Bouwens et al. (2015) is shown at the appropriate 
redshift (z = 12 for the bottom panel), with the shaded region representing the 
l-a uncertainty in the fit parameters. In each panel, solid vertical grey lines 
delineate the detection limits at that redshift for galaxies at the epoch shown 
using (from left to right) the limiting magnitudes of the un-lensed Hubble 
Frontier Fields (mii m = 28.7), the HUDF12 (mii m = 30.1), and the JWST ul¬ 
tradeep campaign (mii m = 31.4). The vertical dashed-grey line corresponds to 
a JWST ultradeep observation with a lensing magnitification of 10. The ma¬ 
genta line in the bottom panel is a toy model of the UV luminosity function 
for this simulation at z = 12, described in Section 4. 
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of the faintest galaxy that could in principle be detected at 
that redshift in an un-lensed Hubble Frontier Field, the Hub¬ 
ble Ultra-Deep Field, and the planned JWST ultradeep field, 
as well as the limiting magnitude for a JWST observation as¬ 
suming a field with a magnification factor of /i = 10. For all 
observational estimates we assume a uniform K-correction of 
-2 magnitudes. In all simulations the UV luminosity function 
matches observations reasonably well at the bright end, but 
there is a flattening of the luminosity function in dim galax¬ 
ies - galaxies that currently cannot be observed, but which in 
the future may be probed by JWST. We also see that the UV 
luminosity function has a termination point, with no galaxies 
dimmer than Muv — -2. This behavior is captured reason¬ 
ably well by a simple toy model, shown in the bottom panel 
and described in Section 4. 

Figure 2 displays the relationship between a given galaxy’s 
stellar mass and the halo’s virial mass or ultraviolet luminos¬ 
ity, combined at their final redshift. In the left panel, a 2D 
histogram of stellar mass versus halo virial mass, lines of con¬ 
stant stellar mass fraction (defined as f* = M*/M v i r ), ranging 
from f* = 10 _1 to 10 -4 , are plotted. The purple solid line and 
accompanying equation is a fit to the median stellar fraction 
as a function of virial mass, 


/» = 1.26 x HP 3 


M, 


V io 8 m 0 


0.74 


( 1 ) 


which is made by fitting a line to the distribution of log stellar 
mass and log virial mass for halos with 10 7 < M v i r /M 0 < 
10 8 - 5 . At virial masses larger than 10 85 M 0 , there is some 
deviation from this relationship, with the stellar fraction in the 
most massive halos being lower than this power-law fit would 
predict. The red-dashed line indicates the stellar fraction as 
a function of virial mass predicted by Behroozi et al. (2013) 
but extrapolated to z = 15 - a fit that agrees with the simulated 
galaxies to a remarkable degree, including the flattening of 
the stellar fraction at high mass. The right panel shows a 2D 
histogram of halo virial mass versus absolute UV magnitude 
for the same galaxies, and indicates a positive (albeit noisy) 
correlation between halo mass and star formation rate. The 
quantities in both panels are insensitive to redshift during the 
epochs considered by this work. 

Figure 3 shows the fraction of halos containing metal- 
enriched stellar populations as a function of halo virial mass 
for all three simulations at their final redshift. All halos 
with M v i r > 5 x 10 7 M 0 contain stars, and all halos with 
M v i r >2x 10 8 M 0 have formed stars in the last 20 Myr (and 
thus are emitting significant ultraviolet light). No halos with 
virial masses below ~ 7 x 10 6 M 0 contain stars. This fig¬ 
ure also contains a fitting function for the fraction of halos 
that have active stars (modeled on the form used by Okamoto 
et al. 2008): 


/occ(M) = 





( 2 ) 


with M c = 6.0 x 10 7 M 0 and a = 1.5. 

Figure 4 shows projections of matter overdensity and elec¬ 
tron fraction in the refined simulation subvolumes at their final 
redshifts. In each panel of the top row, the black bar indicates 
the physical size subtended by a field of view 0.5' on a side 
at that redshift. The projection area is slightly smaller than 
the 2.2' x 2.2' JWST NIRCam field of view. In both rows, all 


galaxies with Mi6oo < -10 are shown, with glyphs showing 
galaxies which could be seen in the 2012 Hubble Ultra Deep 
Field (red triangles) (Koekemoer et al. 2013), the un-lensed 
(blue squares) and lensed (green circles; /i = 10) JWST ultra¬ 
deep field (Gardner et al. 2006), and dimmer galaxies (black 
dots). The spatial distribution of galaxies in each field demon¬ 
strates cosmic variance. The rare peak field shows many clus¬ 
tered galaxies above the lensed JWST limit, whereas the nor¬ 
mal field has very few galaxies outside of a clustered region. 
The void region is bounded on the right by an overdense sheet 
containing several galaxies above the JWST unlensed sensi¬ 
tivity limit. In the void, there is little collapsed structure that 
can sustain efficient star formation, leading to a lack of galax¬ 
ies with Mi goo < -10. The projected electron fraction in the 
lower row depicts the H II regions with radii up to ^1 co¬ 
moving Mpc associated with these galaxies. Scattered around 
the volume are relic H II regions created by Population III 
stars that form in substantially smaller halos (M v i r 10 6 M 0 ) 
than the galaxies shown in the Figure, marginally adding to 
the overall ionized fraction. The overall ionized morphology 
is consistent with an inside-out scenario, previously noted by 
several groups (e.g. Gnedin 2000; Trac et al. 2008; Finlator 
et al. 2009; So et al. 2014) during the initial phases of reion¬ 
ization. This is most apparent in the void simulation, because 
the ionization front from the overdense sheet has not yet prop¬ 
agated into that region and there are no void galaxies to reion¬ 
ize their surroundings. 

4. SUMMARY AND DISCUSSION 
The key results presented in this paper are as follows: 

1. The ultraviolet luminosity function of our simulated 
galaxies is consistent with observations of high-z 
galaxy populations at the bright end of the luminosity 
function (Misoo ^ -17), but at lower luminosities is es¬ 
sentially flat rather than rising steeply as inferred by a 
Schechter function fit to observational data, and has a 
clearly-defined lower limit in UV luminosity. 

2. This behavior of the UV luminosity function is due to 
two factors: (i) the strong dependence of the star for¬ 
mation rate on halo virial mass in our simulated galaxy 
population, with lower-mass halos having systemati¬ 
cally lower star formation rates and thus lower UV lu¬ 
minosities; and; and (ii) the fact that halos with virial 
masses below ^2x 10 8 M 0 do not universally con¬ 
tain stars, with the fraction of halos containing stars 
falling with decreasing virial mass and reaching zero 
at~7xlO 6 M 0 . 

3. The brightest of our simulated galaxies may be visible 
to current and future ultra-deep space-based surveys at 
z ~ 12 (and likely at lower redshifts as well), particu¬ 
larly if lensed regions are chosen for observation. 

The primary result of this paper - the flattening of the mass 
function - is significant for our understanding of the reion¬ 
ization of the universe. Observations of high redshift galax¬ 
ies provide poor constraints on the low-luminosity end of the 
galaxy luminosity function, and thus make it challenging to 
accurately account for the full budget of ionizing photons dur¬ 
ing that epoch. Our work suggests that there are far fewer 
faint galaxies at high redshift than one would infer from fit¬ 
ting a Schechter function to the high-luminosity end, as has 
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Figure 2. Left panel: 2D histogram of galaxy stellar mass versus halo virial mass for all simulations combined at their final redshifts. The histogram is 
logarithmic, and the colors indicate the number of galaxies in each bin (with red/yellow indicating a large number of galaxies and green/blue indicating few 
galaxies). Black-dashed lines indicate constant stellar fractions (f* = M*/M v i r ), ranging by factors of 10 from f* = 10 _1 (top) to 10 -4 (bottom). The purple solid 
line and accompanying equation indicates the fit to the median stellar fraction as a function of virial mass for halos with 10 7 < Mv^/Mq < 10 8 5 . The red-dashed 
line indicates the stellar fraction as a function of virial mass predicted by Behroozi et al. (2013) but extrapolated to z - 15. Right panel: Halo virial mass versus 
absolute UV magnitude for the same set of galaxies. 



Figure 3. Fraction of halos containing stellar populations as a function of 
halo mass for all simulations at their final redshift. Solid lines: fraction of 
halos of that virial mass containing metal-enriched stars of any age. Dashed 
lines: fraction of halos of that virial mass that have experienced star formation 
in the past 20 Myr. Grey solid line: fitting function used in the toy model 
shown in Figure 1 . 

been done by, e.g., Bouwens et al. (2015). Taken at face 
value, a smaller number of galaxies would result in an un¬ 
der-production of photons relative to what is needed to reion¬ 
ize the universe. However, the Planck 2015 results (Planck 
Collaboration et al. 2015), particularly the Thomson scatter¬ 
ing optical depth r = 0.066 ±0.012 corresponding to an in¬ 
stantaneous reionization redshift Zr = 8.8t};?, relieve the pre¬ 
vious tension between reionization constraints from polar¬ 
ization measurements of the CMB and galaxy observations 
(Robertson et al. 2015). 

There are many complicating factors when computing a 
reionization history from a galaxy luminosity function - in 
particular, the UV photon escape fraction and its relationship 
to the halo virial mass (see, e.g., Wise et al. 2014; Fontanot 
et al. 2014) - and we defer a detailed analysis to future work 
(Xu et al. 2015, in prep.). Recent work has shown that low- 
mass galaxies (M v i r < 10 8 M 0 ; M* < 10 5 M 0 ) prior to reion¬ 
ization have typical escape fractions in the range 5-50% with 
the upper range occurring when supernova blastwaves and 


ionization fronts create low-density channels that allow ion¬ 
izing radiation to freely escape into the neutral intergalactic 
medium (Paardekooper et al. 2013, 2015; Wise et al. 2014). In 
more massive halos, the simulations of Kimm & Cen (2014) 
showed that the time-averaged UV escape fraction for a par¬ 
ticular galaxy is only ^10% (although Ma et al. (2015) sug¬ 
gest that this fraction may be considerably lower over halos 
with a mass range of 10 9 - 10 11 M 0 , with a time-averaged 
value of around ~ 5%). This suggests that low-luminosity 
galaxies may provide a non-negligible amount of the ionizing 
photon budget, given their higher escape fractions and num¬ 
ber densities. 

The flattening of the high redshift galaxy luminosity func¬ 
tion at the faint end is a robust result 2 in our simulations (and 
has been suggested by other, albeit less physics-rich, simu¬ 
lations - see Jaacks et al. 2013). What is less robust, how¬ 
ever, is what can be inferred from the overall normalization of 
the luminosity function from our calculations. Given JWST’s 
small field of view, cosmic variance will present a significant 
challenge (e.g., Trenti & Stiavelli 2008). We have sampled 
regions of varying mean density, which gives some sense of 
cosmic variance. However, a precise discussion of this subject 
is difficult because of the small cosmological volumes sam¬ 
pled, as well as the somewhat arbitrary nature of our choice 
of resimulated regions. 

We note that the strongly evolving stellar mass fraction as 
a function of halo mass has an intriguing similarity with ob¬ 
servations of Local Group dwarf galaxies (e.g., McConnachie 
2012), where there is a strong positive correlation between 
a given galaxy’s mass-to-light ratio and its bolometric lu¬ 
minosity. The smallest local dwarf galaxies are strongly 
dark matter-dominated, having similar stellar fractions to our 
lower-mass galaxies, and the trends are similar as well. While 
it is hard to make precise comparisons between Local Group 
observations and high-redshift galaxy simulations, the sim¬ 
ilarity in trends suggests that there may be universal, halo 
mass-dependent behaviors in galaxy formation. 

2 See Supplementary Online Materials. 
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Figure 4. Projections through the refined subvolumes of the Rare peak. Normal, and Void simulations at their final redshifts. Top row: matter overdensity. 
Bottom row: electron fraction weighted by gas overdensity. In each panel in the top row, the black bar indicates the physical size subtended by a field of view 
0.5 X at that redshift, and the projection area is slightly smaller than the 2.2' x 2.2' JWST NIRCam field of view. In both rows, all galaxies with M^oo < -10 
are shown, with those that in principle could be observed by the HUDF12 shown as red triangles, by the unlensed JWST ultradeep field as blue squares, and 
considering a magnification factor of p = 10 in the JWST ultradeep field as green circles, with glyph size proportional to luminosity. Note that any galaxy 
observable in the HUDF could also be seen in the JWST ultradeep field. Galaxies with Migoo < -10 and unobservable by any of these campaigns are shown as 
black dots. 


Figure 3 indicates something intriguing - namely, that there 
is some physical effect that quenches star formation in low 
mass halos in the early universe. The cause is not reionization 
- Figure 4 demonstrates that the universe is patchily reionized 
at this epoch - but rather F^-photodissociating Lyman-Werner 
radiation (e.g., O’Shea & Norman 2008). The universe is 
transparent to this radiation, and thus star formation anywhere 
in the volume can suppress molecular hydrogen elsewhere, 
with the level of suppression growing with the comoving star 
formation rate density. The consequence of this suppression 
is that low-mass halos of primordial composition, which cool 
primarily by H 2 line emission, must grow to be more mas¬ 
sive and hotter before the gas in the halo core can cool effi¬ 
ciently, which suppresses Population III star formation (and, 
thus, metal-enriched star formation). Once halos reach this 
threshold mass (which depends on the strengh of the Lyman- 
Werner background; see O’Shea & Norman 2008), stars can 
then form. We note that at later times, when the universe is 
transparent to hydrogen-ionizing photons, a similar but more 
pronounced effect must exist due to this photon population. 

Figures 2 and 3 strongly suggest that the flattening of the 
galaxy luminosity function is due to a combination of de¬ 
creased star formation rate at low halo virial masses and the 
lack of stars in halos below ~ 10 8 M 0 . A simple toy model 3 

3 See Supplementary Online Materials for a detailed description. 


using the analytic halo mass function convolved with a con¬ 
stant specific star formation rate and the halo stellar occupa¬ 
tion fraction from Figure 3 is shown as the magenta line in 
the bottom panel of Figure 1, and matches the properties of 
the luminosity function of our simulated galaxies quite well, 
accurately capturing the deviation from a Schechter function, 
the overall flattening at intermediate luminosity, and the ter¬ 
mination at the low-luminosity end. 

The dependence of stellar fraction on halo mass is well- 
known observationally (e.g., Behroozi et al. 2013), and is also 
seen in simulations of more massive galaxies at lower red- 
shifts (e.g., Hopkins et al. 2014). We speculate that this be¬ 
havior is due to a combination of cooling and stellar feed¬ 
back (as has been discussed in a different context by Voit 
et al. 2014). Gas in low-mass halos cools inefficiently due to 
the low virial temperature, and stellar feedack can efficiently 
remove metal-enriched gas from these low-mass halos (see, 
e.g., Chen et al. 2014), which further decreases the efficiency 
with which gas can cool and form stars. The suppression of 
star formation in halos below ~ 10 8 M Q (and complete ab¬ 
sence below ~ 10 7 M 0 ) is likely due to inefficient cooling, 
destruction of H 2 via Lyman Werner radiation from neighbor¬ 
ing halos, and the resulting Jeans screening (e.g., Gnedin & 
Hui 1998; Wise & Abel 2008; Wise et al. 2014). This sug¬ 
gests that there is an effective lower limit on the masses of 
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high-redshift galaxies, and thus a definitive lower end to the 
UV luminosity function. We will explore this in future work. 
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APPENDIX 

SUPPLEMENTARY ONLINE MATERIAL 

A Toy UV Luminosity Function Model 

The simple luminosity function model shown by the magenta line in the bottom panel of Figure 1 was created using the 
following procedure: 

1. Calculate the comoving number density of cosmological halos as a function of halo virial mass m, at z = 12 for the 

WMAP7 best-fit cosmology, using the fitting function of Warren et al. (2006). 

2. For each halo mass bin, calculate the recent star formation history of the halo assuming a constant specific star formation 
rate of 10 -8 yr -1 over the last 100 Myr. This results in a constant star formation rate, under the assumption of constant halo 
mass. 

3. Given this star formation history and the assumption of a simple Salpeter IMF with parameters identical to those in the z = 0 
Milky Way (consistent with our simulation treatment of metal-enriched stars; see Section 2), calculate the IMF-integrated 
bolometric luminosity of the stars that have not yet undergone supernovae at the epoch under consideration. Given that 
this luminosity is dominated by the most massive and hottest stars (since luminosity scales approximately as stellar mass 
to the fourth power), we make the assumption that the UV luminosity of the stellar population in each halo is roughly its 
bolometric luminosity. This results in a comoving number density of halos as a function of UV magnitude, d y ^ M , where 
M represents absolute UV magnitude. 

4. Finally, we multiply the resulting quantity by the fraction of halos that have experienced recent star formation as a function 
of halo mass that is shown in Figure 3, using the crude approximation shown by the grey line. This fit captures the essence 
of the figure - namely, universal recent star formation above M v i r >2x 10 8 M 0 , no recent star formation below rsj 10 7 Mq, 
and a steady decrease in star formation with decreasing halo mass between those endpoints. 




Probing the High-z UV Luminosity Function 


7 


This toy model captures several of the essential features of the luminosity function shown by our simulations: it has a significant 
slope at the high-luminosity end, flattens out at Migoo > -12, and goes entirely to zero at M^oo — -2. Furthermore, the properties 
of this toy model are robust to several of the simplifying assumptions made here. The most critical of these are as follows: 
Specific star formation rate: The specific star formation rate (SSFR) evolves as both a function of redshift and halo mass. This 
evolution is most pronounced at low redshift (z < 2) - at higher redshifts, such as those considered in this work, the SSFR evolves 
very slowly with redshift and is relatively insensitive to halo mass (see, e.g., Appendix F of Behroozi et al. 2013), making the 
approximation of a constant SSFR sensible. We note that at the redshifts and halo masses of interest to the work presented here 
this quantity is extremely hard to constrain observationally; however, extrapolating the observed trends to z > 8 suggests little 
evolution. Deviation from a constant SSFR will modify the predicted UV luminosity function in proportion to the deviation. 
Halo formation history: Cosmological halos grow over time, and do so particularly quickly at the redshifts considered in this 
work. Assuming that the specific star formation rate is directly proportional to halo mass, this implies that the total (as opposed to 
specific) star formation rate in a given halo is increasing rapidly as well. The period of 100 Myr over which the star formation rate 
is assumed to be constant is a significant fraction of the Hubble time at high redshift (~ 250 Myr at z = 16, ~ 380 Myr at z = 12), 
suggesting that this assumption may be suspect. However, the UV luminosity is dominated by massive stars having lifetimes of 
millions of years, and thus the production of UV photons is dominated almost entirely by the very recent star formation history 
of the halo. As a result, the UV luminosity is quite robust to the exact details of the halo’s formation history, with deviations from 
constant star formation rate contributing mildly to variation in magnitude. 

Comparison with Higher Resolution Simulations 

The simulated behavior of the UV luminosity function at the faint-end can be suspectible to limited numerical resolution. To 
check for any spurious effects, we compare our results to our previous work (Wise et al. 2012b, a, 2014) that focused on 32 
first-generation galaxies in a small comoving volume of (1 Mpc) 3 with a DM particle mass resolution of 1840 M 0 , 15 times 
smaller than the dark matter particles in the simulations presented in this paper. The major shortcoming of our previous work 
was the limited galaxy sample size, which was rectified in the Renaissance Simulations. These two simulation suites share 
the same subgrid physics models — chemical reaction networks, metal-free and metal-enriched star formation, and radiative 
and supernova feedback. Although the Renaissance Simulations cannot capture Population III star formation in the smallest 
minihalos (^2x 10 5 M 0 ), properties of metal-enriched star formation, in particular the stellar mass - halo mass relation, agree 
very well between the two simulation suites (see Fig. 4 in Chen et al. 2014). There we showed that this relationship in both 
simulations were within 1-cr of each other in halos with M v i r > 10 7 M 0 , which are the smallest halos considered in our analysis. 
Furthermore, Wise et al. (2014) found that the UV luminosity function was flat at magnitudes M^oo ^ —12, consistent with our 
findings in the Renaissance Simulations. Thus, we are confident that the flattening of the luminosity function is not a result of 
limited numerical resolution but is caused by the suppression of star formation from radiative and supernova feedback effects in 
halos with masses M v i r <2x 10 8 M 0 . 


